Scaling Laws and Temperature Profiles for Solar and Stellar 
Coronal Loops with Non-uniform Heating 
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ABSTRACT 

The bulk of solar coronal radiative loss consists of soft X-ray emission from quasi-static loops 
at the cores of Active Regions. In order to develop diagnostics for determining the heating 
mechanism of these loops from observations by coronal imaging instruments, I have developed 
analytical solutions for the temperature structure and scaling laws of loop strands for a wide range 
of heating functions, including footpoint heating, uniform heating, and heating concentrated at 
the loop apex. Key results are that the temperature profile depends only weakly on the heating 
distribution - not sufficiently to be of significant diagnostic value - and that the scaling laws 
survive for this wide range of heating distributions, but with the constant of proportionality in 
the RTV scaling law {P^L ~ T^^^) depending on the specific heating function. Furthermore, 
quasi-static analytical solutions do not exist for an excessive concentration of heating near the 
loop footpoints, a result in agreement with recent numerical simulations. It is demonstrated 
that a generalization of the solutions to the case of a strand with a variable diameter leads to 
only relatively small correction factors in the scaling laws and temperature profiles for constant 
diameter loop strands. A quintet of leading theoretical coronal heating mechanisms is shown to 
be captured by the formalism of this paper, and the differences in thermal structure between 
them may be verified through observations. Preliminary results from full numerical simulations 
demonstrate that, despite the simplifying assumptions, the analytical solutions from this paper 
are stable and accurate. 

Subject headings: Sun: chromosphere — Sun: corona — hydrodynamics 



1. Introduction 

The emission from the solar corona in the 
Extreme Ultra- Violet (EUV) and in X-rays has 
been observed since the advent of space-borne 
telescopes. The mere existence of that emission 
demonstrates that the solar corona has a very 
high temperature and therefore is not in radiative 
equilibrium with the underlying photosphere. The 
solar corona requires energy deposition of some 
sort to maintain its high temperature, and it is 
now widely assumed that this coronal heating is 
magnetic in nature, either through the dissipation 
of electrical currents or of incoming MHD waves. 
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The S-054 X-ray telescope (Vaiana et al. 1973) 
on the Skylab mission in particular has contributed 
greatly to elucidating the nature of the coronal 
high temperature emission. It became clear that 
the emission from the corona is highly structured 
in so-called coronal loops; elongated bent high as- 
pect ratio cylinders that appear to be rooted in 
the solar chromosphere at both ends. The Skylab 
observations further showed that the footpoints 
of the X-ray emitting coronal loops are predomi- 
nantly located at photospheric plage regions or the 
penumbrae - but not at the umbrae - of Sunspots. 
This demonstrates a clear connection with en- 
hanced photospheric magnetic field, but curiously 
not with the strongest field. After Skylab soft X- 
ray images from the solar corona have been exten- 



sively studied by the Soft X-ray Telescope (SXT) 
onboard Yohkoh (Tsuneta et al. 1991) and cur- 
rently by the X-ray Telescope (XRT, Golub et al. 
2007) onboard the follow-up Japanese solar mis- 
sion, Hinode. The latter has the unprecedented 
resolution in X-rays of two arcseconds, a factor of 
2.5 or more better than its predecessors. 

While the coronal soft X-ray emission origi- 
nates from loop plasma at typical temperatures 
of several (2-5) MK, solar EUV emission comes 
from plasma with a temperature of 1-2 MK. The 
EUV full disk corona has been observed for an 
entire solar cycle now by the Extreme Ultravio- 
let Telescope (EIT, Delaboudiniere et al. 1995) 
onboard the Solar and Heliospheric Observatory 
{SoHO), and since 1998 with the unparalleled res- 
olution of one arcsecond by the Transition Re- 
gion and Coronal Explorer ( TRA CE, Handy et 
al. 1999). The TRACE observations in particu- 
lar show that the EUV corona at ~ 1 MK (in the 
171 A passband) is clearly structured, with much 
more pronounced loops than the hotter corona, 
as observed in narrow-band images, for example, 
in the TRACE 285 A passband. The same ap- 
pears, interestingly enough, not to be so in broad- 
band coronal soft X-ray images, such as those from 
XRT. I point out here that the TRACE images 
are usually shown on a linear intensity scale after 
subtraction of the dark current, while soft X-ray 
images, such as those from SXT and XRT, usu- 
ally show the square root of the image intensity 
or the natural logarithm, also after subtraction 
of a pedestal. The latter rendition is intended 
to bring out the fainter structures more clearly, 
motivated by the fact that in soft X-rays coro- 
nal emission tends to be dominated by a few very 
bright loops, i.e., the dynamic range in intensity 
of the emitting structures is much larger and dif- 
ferently distributed than in the EUV. This man- 
ner of image rendering also has the unintended 
effect of enhancing the background in soft X-ray 
images in comparison with EUV images, adding 
to the appearance of more "fuzzy" coronal images 
that is not real. In fact, recent work by Cirtain et 
al. (2006) demonstrates that in typical Active Re- 
gions more than three quarters of the emission in 
the EUV originates from the unresolved TRACE 
background (i.e., non-loop emission), while inspec- 
tion of the Hinode soft X-ray images yields a frac- 
tion of a fifth at the very most (Weber 2008, pri- 



vate communication). 

The relevance of the above is that the soft 
X-ray and EUV emission from the solar corona 
is completely dominated by that of solar Active 
Regions; Aschwanden (2001) and Aschwanden et 
al. (2007) estimate that the Active Regions ac- 
count for 82.4 % of the heating requirement of 
the corona. Therefore, in order to understand the 
origin of solar coronal heating one needs to con- 
sider the heating of Active Regions (ARs), and in 
particular the heating of soft X-ray emitting AR 
loops, which dominate the energy output. 

Typically one will find soft X-ray loops in the 
center of ARs, and the somewhat cooler and of- 
ten longer EUV loops in the periphery. In addi- 
tion there is substantial EUV emission from the 
so-called "moss", located at and just above the 
footpoints of the hotter soft X-ray loops. This is 
illustrated by Figure [U showing two simultaneous 
and co-aligned coronal images, one in the EUV 
taken by TRA CE and one in soft X-rays taken by 
XRT on Hinode. 

The hotter soft X-ray loops typically appear 
relatively steady over time-scales that are longer 
than their radiative and conductive cooling times, 
while the EUV loops seem more dynamic, often 
exhibiting a cooling collapse as well as jets and 
siphon flows of chromospheric material entering 
these loops. The absence of significant EUV emis- 
sion, except for the moss, in the soft X-ray loops 
at the centers of ARs led Antiochos et al. (2003) 
to the conclusion that these loops must be either 
steadily heated, or intermittently - like, say, in 
nano-flares - with a repeat time that is shorter 
than the loop cooling time: if that weren't the 
case these loops should be observed emitting in 
EUV bands as they cool down. Recent work by 
Winebarger, Warren, & Falconer (2008) confirms 
this result. Antiochos et al. quantify their result 
with numerical simulations of loop hydrodynam- 
ics using the sophisticated NRL loop code. They 
further demonstrate that since the observed moss 
emission from the footpoints of the soft X-ray 
emitting loops varies little, of the order of 10% 
over hours, the coronal loop heating mechanism 
must be rather steady in time. This justifies the 
quasi-static modeling of the thermal structure of 
the AR soft X-ray loops, as will be done in the 
present paper. Quasi-static is defined here as hav- 
ing the constraints on the heating rate mentioned 
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Fig. 1. — Simultaneous (within one sec) images of an Active Region (AR) with the XRT Thin Be filter, and 
the TRACE 171 A filter. The arrows, at identical locations in both images, point to, from left to right and 
top to bottom: a) a fan of EUV loops in the periphery of the AR, b) hot Soft X-ray loops in the core of the 
AR, with EUV moss at their footpoints, c) a second fan of EUV loops from near the umbra of a sunspot, 
with no corresponding soft X-ray emission, d) diffuse hot and very long Soft X-ray loops surrounding the 
core, with no corresponding EUV emission, and e) short bright soft X-ray loops with corresponding EUV 
moss. Image courtesy of Dr. Leon Golub. 



above, and having plasma-flows in the loops that 
are substantially subsonic (say, less than 1/3 of 
the speed of sound) . The analysis of recent obser- 
vations with XRT and EIT on Hinode by Warren 
et al. (2007) confirms this picture by showing a 
steady output for the bulk of the emission from the 
AR analyzed, supplemented by transient brighten- 
ings in a limited number of loops. 

It was mentioned above that the EUV loops 
in the periphery of ARs appear much more dy- 
namic. Yet Aschwanden, Schrijver, & Alexander 
(2001) find from the analysis of some 40 loops ob- 
served by TRACE that about 30% of these EUV 
loops appear in quasi-steady hydrostatic equilib- 
rium. For the remainder that is dynamic one can 
show that a quasi-static model is a good zeroth 



order approximation for the thermal structure as 
long as any plasma motions are substantially sub- 
sonic, as defined above. The reason for that is that 
the kinetic terms in the energy equations scale as 
the square of the ratio of the plasma velocity to 
the sound speed. Therefore, also for the more dy- 
namic EUV AR loops quasi-static modeling is a 
useful first approach, certainly if one can obtain 
general analytical results as I will do in this pa- 
per. 

Having argued that coronal heating is concen- 
trated in soft X-ray emitting AR loops, and that 
for these loops quasi-static modeling of the ther- 
mal structure is an acceptable approach, I will 
derive in the present paper full analytical solu- 
tions for the thermal structure of coronal loops 



based upon these premises. Analytical solutions 
and approximations have been derived by many 
since the pioneering papers by Rosner, Tucker, 
& Vaiana (1978b) and Craig, McClymont, & Un- 
derwood (1978) but full analytical solutions have 
been limited to the case of uniform quasi-static 
heating of loops. Numerical simulations, on the 
other hand, have investigated non-uniform heat- 
ing, in particular cases where heating is concen- 
trated near the loop footpoints (e.g., Serio et al. 
1981; Aschwanden et al. 2001; Miiller, Hansteen, 
& Peter 2003; Miiller, Peter, & Hansteen 2004; 
Miiller et al. 2005). Here I will derive analyti- 
cal solutions for the loop temperature profile and 
scaling laws for non-uniform heating functions, in- 
cluding cases in which the heating is concentrated 
near the footpoints and around the looptop. These 
solutions produce generalizations of the now 30 
year old RTV scaling laws for uniform heating 
that have some surprising properties that, as I will 
show, may be used as a diagnostic for determin- 
ing the heating profile in the loops. On the other 
hand, it will be demonstrated that the tempera- 
ture profile of coronal loops is rather insensitive 
to the distribution of the coronal heating, so that 
even very precise measurements of the tempera- 
ture profile along loops would not be very useful 
for determining the heating distribution. 

In modeling coronal loops it is usually assumed 
that the loop has a uniform cross-section along its 
length. This is what observations, both in EUV 
and soft X-rays, indicate to be the case (Klimchuk 
2000, Watko & Khmchuk 2000). Visual inspec- 
tion of Figure [1] seems to confirm this. However, 
a recent paper by DeForest (2007) has cast doubt 
upon this notion. To understand the argument I 
first need to introduce the concept of "strands", 
a term that will be used throughout this paper. 
Conduction across the magnetic field is very in- 
effective in coronal plasma. The basic reason is 
that the mean free path for perpendicular particle 
motions is limited to the ion gyro radius, which is 
only of the order of centimeters, while the mean 
free path along the magnetic field can be hundreds 
of kilometers or more for thermal electrons. Hence 
conductive energy transport along the magnetic 
field is far more effective than across the field. 
That results in the situation, already recognized in 
Rosner et al. (1978b) and used ever since in loop 
modeling, that the solar corona is divided in inde- 



pendent mini-atmospheres, or "strands", stretch- 
ing out along magnetic field lines from one chro- 
mospheric footpoint to the other, but with very 
tiny cross-sections, possibly down to an ion gyro- 
radius. These strands are then modeled as one- 
dimensional structures, with the single coordinate 
along a field line. 

The loops observed by solar imagers, from the 
Skylah telescopes to TRACE and XRT on Hin- 
ode, although having aspect ratios of 10% and 
smaller, still have radii of several thousands of 
km, down to perhaps 1000 km for the thinnest 
TRACE loops. Therefore these observed loops 
may consist of many, up to thousands, of indi- 
vidual "strands". Note that it is interesting, and 
perhaps a clue to the physical mechanism of the 
heating mechanism, that these "stands" tend to 
cluster together in the much wider loops that are 
observed - there is no a priori reason strands could 
not be distributed evenly over ARs, thereby cre- 
ating a much smoother appearance. The many 
strands within an observed loop can all have dif- 
ferent heating rates and temperature structures. 
On the other hand, it is possible that the observed 
loops consist of many strands that are heated in 
an identical manner and therefore all have the 
same temperature structure. It was pointed out by 
Martens, Cirtain, & Schmelz (2002) that if the for- 
mer is true one would expect that pixels from these 
loops have plasma from different strands at multi- 
ple temperatures contributing to the emission, and 
therefore that the Differential Emission Measures 
(DEMs, Pottasch 1963) for these pixels would be 
broad. Martens et al. then showed that for broad 
loop DEMs the customary method of determining 
temperatures from filter ratios (for narrow-band 
images) and line-ratios (for spectrometers) leads 
to misleading results. They demonstrated that for 
broad DEMs the ratio for one pair of lines or pass- 
bands in general produces a different temperature 
than those from another pair, exactly what has 
been observed (e.g., Schmelz et al. 2001; Reale 
et al. 2007). Weber et al. (2005) carried this ap- 
proach further through simulations and found just 
how narrow a DEM has to be for the filter or line- 
ratio method to yield a single, correct, filter-ratio 
temperature. In general these results imply one 
cannot assign a single temperature to each pixel 
within a multi-strand multi-temperature loop. 

The answer to the question whether observed 



loops are multi-thermal or isothermal, obviously 
depends on the resolution of the telescope that 
is used: once we get down in resolution to the 
size of the elementary strands, the answer must 
be isothermal. However, so-far the answer for Soft 
X-ray loops observed by the currently highest reso- 
lution soft X-ray telescope, XRT, is multi-thermal 
(Reale et al. 2007), and therefore the observations 
can only be matched by models that contain a 
multitude of elementary strands. This leaves the 
modeler with a large number of parameters to ad- 
just (number of strands and heating rate and pro- 
file in each) which implies analytical solutions for 
each individual strand - as will be derived in this 
paper ~ are very useful, since it makes it much 
less computationally demanding to explore the pa- 
rameter space of the models in trying to obtain a 
forward fit to the data. For EUV loops there is a 
heated debate on whether the thinnest loops ob- 
served by TRACE are isothermal (e.g., Aschwan- 
den & Nightingale 2005) or multi-thermal (e.g.. 
Martens et al. 2002; Schmelz & Martens 2006; 
Schmelz, Kashyap, & Weber 2007), an issue that 
I will not further address in this paper. The ob- 
servations by the Atmospheric Imaging Assembly 
(AIA) onboard SDO, to be launched by the end of 
2008, might resolve this question for the TRACE 
loops, since AIA has the same spatial resolution 
as TRACE, but a much better spectral resolution 
due to its larger number of EUV passbands. Re- 
gardless of the outcome of the debate, it would 
be a great step forward when the first isother- 
mal strands are detected, because that puts strong 
constraints on the physics of the heating mecha- 
nism. 

As an aside, from the considerations above it 
also follows that the loops observed in the EUV 
and soft X-rays also outline the magnetic field, 
and therefore are a great diagnostic for the struc- 
ture and topology of the coronal magnetic field, as 
the high resolution TRACE images so beautifully 
demonstrate. 

Now let's return to the issue of the cross-section 
of loops and strands. In order to model a strand 
one has to make an assumption about the vari- 
ation of its cross-section, i.e., is it constant or 
does the cross-section increase from footpoint to 
top? As was stated above, observed loops tend 
to show a constant cross-section, but that is not 
necessarily true for the strands within them, as 



pointed out by DeForest (2007) . Extrapolations of 
the coronal magnetic field from the photospheric 
boundary nearly always yield flux tubes expand- 
ing with height, so why is it that observed loops 
do not follow that trend? DeForest shows a clear 
example of a polar plume observed by FIT (his 
Fig. 3) that appears to have a constant diame- 
ter, but the same plume observed by TRACE at a 
five times higher resolution, consists of a number 
of smaller strands, while the whole of the plume 
clearly shows an increasing cross-section from the 
bottom up. The analysis by DeForest demon- 
strates that near the limit of its resolution any 
telescope will show strands that are expanding as 
having a constant diameter, because of the smear- 
ing inherent in the Point Spread Function of the 
telescope, because of photon noise, and because 
of the background and foreground emission along 
the line-of-sight. DeForest then shows that if these 
strands are in fact expanding this can explain 
why the observed emission measure at the top of 
loops in TRACE is typically larger than would 
be expected from modeling constant cross-section 
strands - simply because there is more emitting 
volume near the top. I should point out here that 
others (e.g., Lopez- Fuentes, Demoulin, & Klim- 
chuk 2008) strongly disagree with DeForest 's re- 
sults, arguing, based upon their simulations, that 
even for loops that are close to the telescope's res- 
olution the diameter can be correctly measured. 
However, the authors offer no suggestion how to 
reconcile the observations of constant loop diam- 
eter with magnetic filed extrapolations almost in- 
variably showing the opposite. (To be fair, that is 
outside the scope of their paper). 

This issue not being resolved yet, I have de- 
veloped analytical solutions for strands both with 
constant and varying diameters. It is shown that 
the initial solution derived for constant cross- 
section strands can be fairly easily generalized to 
a class of non-constant diameter solutions, and, 
most importantly, that the resulting temperature 
profiles and scaling laws only differ by a small 
correction factor from the uniform diameter case. 
This fortunate coincidence makes it possible to 
produce forward folding models for multi-strand 
loops without having to choose between uniform 
or varying strand diameter, while at the same time 
this result allows for an observational verification 
of the degree of uniformity of the strand diame- 



ter, because the emitting volume within any loop 
varies proportional to the change in the strand 
cross-section, as already analyzed by DeForest 
(2007). 

In summary of this introduction: in this paper 
I will develop analytical solutions for the temper- 
ature profile and scaling laws of the strands that 
make up the AR solar coronal loops observed by 
EUV and soft X-ray telescopes. The emission of 
these loops represents the bulk of the solar coro- 
nal emission, and the comparison of their thermal 
structure with the results from analytical and nu- 
merical analyses can therefore contribute to un- 
veiling the detailed physics of the coronal heating 
mechanism. 

In §[2] I will first introduce and validate the ap- 
proximations made in rendering the equations for 
the thermal structure of strands amenable to ana- 
lytical treatment. I will then derive explicitly so- 
lutions that apply to a whole range of non-uniform 
heating functions, from footpoint heating to apex 
heating. I will present both the standard case of 
a uniformly increasing temperature from base to 
apex, as well as solutions with temperature rever- 
sals. 

In § [3] I will present a physical analysis of the 
results, and in particular consider their applica- 
bility as a diagnostic for the heating mechanism 
of strands. To illustrate this I will formulate five 
often considered coronal heating mechanisms in 
terms of the expression used for the heating func- 
tion in this paper, and compare the observable 
predictions for the resulting temperature profiles 
and scaling laws. In § 13.51 1 will generalize the so- 
lutions for the temperature structure and scaling 
laws from the case of a uniform loop cross-section 
to the case of a varying one, in particular the sit- 
uation of a gradually expanding loop diameter in 
going from footpoint to looptop. Finally, in § 13.61 1 
will analyze the solutions exhibiting temperature 
reversals, that, I will argue, apply to prominences. 

In the Discussion section I will conclude by de- 
scribing the heritage of the mathematical descrip- 
tion used in this paper, and comment on the sever- 
ity of the approximations used in the analytical 
treatment as well as on the existence, stability, 
and accuracy of the solutions. 



2. Formalism 

The classical derivation of the RTV scaling laws 
(Rosner et al. 1978b) uses a coronal heating func- 
tion that is a constant in space and time as well 
as constant pressure throughout the loop. The 
latter assumption is reasonable when the height 
of the top of the loop above the chromosphere is 
less than the pressure scale height for the peak 
temperature. Landini & Monsignori-Fossi (1975), 
followed by many others including Martens et al. 
(2000), approximate the radiative loss function for 
constant pressure by a single power-law. The lat- 
ter then explicitly demonstrated that this approx- 
imation is quite accurate for the relevant temper- 
ature range of the corona and transition region. 
Kuin & Martens (1982) showed that this leads to 
particularly simple analytical expressions for the 
temperature profile in an isobaric and uniformly 
heated coronal loop and to exact expressions for 
the two scaling laws. 

I will now derive equally simple results for a 
parametrised non- uniform heating function that, 
depending upon the choice of the parameters, can 
represent heating concentrated at the footpoints 
of the loop, at the apex, and heating that is con- 
stant. I will show that while the form of the first 
RTV scaling law (PqL ~ T^^^, with Pq the loop 
pressure, L, the half-length of the loop, and Tmax 
the maximum temperature at the apex) is pre- 
served for all these heating functions, the value 
of the constant of proportionality depends on the 
spatial distribution of the heating. This constant 
becomes larger as heating becomes more concen- 
trated at the footpoints, so that for the same Tmax 
and L the loop pressure and hence average density 
are larger than for the case of constant heating. 
Therefore these loops are, in fact, not overdense 
with respect to the scaling laws, only overdense 
with respect to the scaling law for uniform heat- 
ing. 

2.1. Uniformly Increasing Temperature 
Solution 

For a static coronal loop under these assump- 
tions the momentum equation yields constant 
pressure and the energy equation reduces to (e.g.. 
Martens et al. 2000) 

^J^oT'/'^) + E, p„2xor-(2+^^ = 0. (1) 



Here z is the coordinate along the axis of the 
strand, measured from the base of the transi- 
tion region, T(z) is the temperature profile, E^ 
the unknown heating function, and xo a-nd 7 
the parametrisation constants in the power-law 
approximation of the radiative loss function for 
constant pressure. Using the Feldman (1992) 
abundances that apply to closed coronal regions. 
Martens et al. (2000) derive that good values are 
7 = 0.5, and xo = 10^^'*^ in cgs units. The 
remaining term in this equation represents con- 
ductive heating or cooling, with kq a constant of 
value 1.1 X 10~^ ergcm^^ K^^ (Chapman & Cowl- 
ing 1939; Spitzer 1962). . 

Rather than assuming uniform heating I now 
use a heating function that has a power-law de- 
pendence upon density and temperature E/j ~ 
pf3 rp{a+f3)^ which, as I will demonstrate below, can 
plausibly represent a number of physical coronal 
heating mechanisms. Using the gas-law to elimi- 
nate the density in favor of the constant pressure 
one finds 

Eh - i/Po^T", (2) 



where H is a constant of proportionality. Clearly 
the double power-law dependence upon density 
and heating collapses to a single power-law for con- 
stant pressure after using the gas law, but this is 
not the case when the pressure scale height be- 
comes comparable to the loop height. Uniform 
heating is represented by the parameter choice 
a=/3=0, while for a uniformly increasing temper- 
ature from footpoint to loop apex, negative a's 
represent footpoint heating, and positive a's apex 
heating. This formalism for non-uniform coro- 
nal heating was introduced by Craig, McClymont, 
& Underwood (1978), and applied by Landini 
& Monsignori-Fossi (1981), TorriccUi-Ciamponi, 
Einaudi, & Chiuderi, and Kuin & Martens (1982). 
Bray et al. (1991) and Kano & Tsuneta (1995, 
1996) used it in an abbreviated form, leaving out 
the pressure term (i.e., (3 = 0), which does not 
make a difference in the mathematical treatment, 
but restricts the physical applicability. 

The energy equation is now cast in a dimen- 
sionless form by introducing the variables 
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as well as the parameters 
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where " represents the second derivative with re- 
spect to x, while ', with obvious meaning, will also 
be used below. 

Equation ([9]) is a second order nonlinear ordi- 
nary differential equation, to be solved in the do- 
main X G [0, 1] and r\ g [0, 1] There are two pa- 
rameters e and ^, that will be shown to be related 
to the scaling laws. With the given form of the 
heating function and the radiative loss function 
this equation can be solved analytically, as will be 
done below. 

Before doing that I will first discuss the bound- 
ary conditions. It is usually assumed, guided by 
observations, that the maximum temperature is 
reached at the top. Hence the boundary condi- 
tions at the loop apex (x=l) become 77(1 )=1, and 
?7'(1)=0. Since equation ([9|) is a second order or- 
dinary differential equation these two boundary 
conditions suffice for determining a solution. How- 
ever, the physical nature of the interface between 
the loop and the underlying chromosphere at the 
footpoints produces two more boundary condi- 
tions. The temperature of the chromosphere is 
of the order of IC* K, while a typical solar coronal 
loop has an apex temperature of the order of 1-5 
MK. Hence the boundary condition at the bottom 
satisfies ?7(0) = (T,;„.„„/T„a.)^/' < 10"^. I will 
show below that for small 77(0) the solutions con- 
verge towards the solution for 77(0) = 0. The same 
is true for the dimensionless conductive flux, t]' . 
The vertical extent of the chromosphere is 2000- 
5000 km, so the average chromospheric temper- 
ature gradient is of the order 3 x lO^^Kcm^^. 
Hence, for a loop length of the order of 10^ cm, 

^'(0) - ""^"^^'T ^ 3 X 10-^ and agam 



the solution can be developed as a series expansion 
in that small parameter. Therefore I can apply the 
footpoint boundary conditions ?7'(0) = 77(0) = 0. 

These extra boundary conditions imply that 
usually no solutions exist to equation ([9]) for arbi- 
trary values of the parameters e and ^. In general 
solutions can be found only for one or more dis- 
crete values of these parameters!^. I will show be- 
low by explicit integration that the extra boundary 
conditions plus the requirement of a single peaked 
temperature profile lead to single unique values 
of e and ^. Writing out these expressions explic- 
itly yields the two scaling laws. Hence, from a 
mathematical point of view, the scaling laws are 
the result of the differential equation for the tem- 
perature structure being overdetermined by two 
excess boundary conditions. Physically the scal- 
ing laws express that for a loop with given length 
and heating profile, the temperature structure, in 
particular Tmax and the pressure Po, are fully de- 
termined, a result that has been confirmed many 
times by numerical simulations (e.g., Serio et al. 
1981; Aschwanden et al. 2001; Miiller et al. 2005). 

The solution of equation ^ is straightforward. 
Multiplication on both sides with 77', and applying 
the boundary conditions at x=0, results in the first 
integral, 

' ^' (10) 
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Applying the remaining boundary conditions at 
x=l produces a generalized scaling law, usually 
called the second scaling law. 



e = 



1 



2a + 7 



/i + 1 3 - 27 ' 



(11) 



which shows that ^ takes on a specific numerical 
value for a given heating mechanism (a) and fit to 
the radiative loss function (7). I will demonstrate 
below that this indeed reproduces the second scal- 
ing law in the more familiar form. Inserting the 
formal boundary conditions ry(0) and ?y'(0) intro- 
duces only vanishing corrections to this scaling law 

e(V^-V(0)^) V'"^^ 77(0)''+i e^'^+i ^^7(0) 



1^+1 



H + l n + l v + l 



(12) 



^This is a typical for similar second order differential equa- 
tions with free parameters, e.g. the Schrodinger equation, 
and others of Sturm-Liouville type, such as the Legendre 
and Bessel equations. 



Applying the boundary conditions at x=l the 
second generalized scaling law for arbitrary condi- 
tions at the foot points of the loop is found 

z/ + 1 1 - 7?(0)"+i + er;'(0)2(/i + l)/2 



e 



M 



1 



(13) 



1 - '/o+^ 

For 7=0.5, ^i =-5/7, and for a > -7/2, t/+l > 0. 
Hence there are no singularities in equation (fT3|) 
under these conditions, and the correction terms 
are vanishingly small for realistic chromosphcric 
boundary conditions. 

What remains is the first order nonlinear ordi- 
nary differential equation (jlOp . It can be expressed 
as 



_ 2(r;^+i-r,-+i) ,, 
^ ^ e(u-fl) ^ ' 



(14) 



where the positive root has been taken conform 
with a monotonic rise in 77 (end hence tempera- 
ture) from the bottom to the top of the loop. This 
equation can be directly integrated and yields 

<■' dy 



x= Ka. + 1)/2]1/2 



(yM+i_y-+i)i/2- 



(15) 



This expression can be simplified by the transfor- 
mation 

u = 7]''-^ (16) 

with a similar substitution of y for z in the inte- 
grand. The integral now reduces to 



e(a. + 1)/2]V2 



^dz 



(17) 



v-^i 7o (i-z)i/2' 

where the parameter A is defined as 

A = ^-^i^^- (18) 

The integrand in equation (|17p represents the def- 
inition of the well-known incomplete /3-function 
(e.g., Zwillinger,1996, p. 497), i.e., 

z^dz 



{i-zyi^ 



/5(u;A + l,l/2), 



(19) 



where Re(A) > -1, so that the integral does not 
diverge. In the appendix it will be shown that the 
effect of a non-zero chromosphcric temperature at 
the lower boundary is very minor. Restoring the 
original parameters a and 7, one finds 

A = m^^^:^^ (20) 



2(2 -f 7 + a) ' 
e(l-2(2 + 7)/7) 1/2 7/3(u;A-H,l/2) 
^ 2 ^ 2{2 + -i + a) ' 



(21) 



Thus the condition A > -1 implies a > -(2+7), i.e., 
a > -2.5 for 7=0.5. 

The boundary conditions at the bottom of the 
loop are satisfied by construction, and applying 
the boundary conditions at the top of the loop 
yields 



1 



e(l-2(2 + 7)/7) 1;, 7i?(A + 1,1/2) 
^ 2 ^ 2(2 + 7 + a) 



(22) 



where B is the complete /?- function (i.e., the inte- 
gral carried out up to u=l). This condition sets 
the value of the parameter e, which, as I shall show 
below, defines the first (RTV) scaling law. 



(l-2(2 + 7)/7) 



2(2 + 7 + a) 
[ 7 J 



1 



S2(A-M,l/2) 



Eliminating e from equation (PT|) yields 

X = /3r(u;A + l,l/2), 



(23) 



(24) 



where (3r is the regularized incomplete /3-function 
(i.e., divided by the corresponding complete 
value). 

Inverting this equation and restoring the origi- 
nal dimensional variables produces an explicit ex- 
pression for all temperature profiles defined by 
valid choices of the parameters a and (3 in the 
heating function Eh = HP^T°', 

T{z) = T,„,,[/3-i(s/L;A + l,l/2)]3TW^, (25) 

where (3~ ^ is the inverse of the regularized incom- 
plete /3- function. Note that the parameters H, Pq, 
and P from the heating function do not appear in 
this solution. I will show below that H and Pq do 
show up in the scaling law relations for Tj^ax i ^^nd 
hence that they do influence the solution. 

2.2. Solutions with Temperature Rever- 
sals 

When the condition of a uniform increase 
in temperature is dropped but the condition of 
77' = at the boundaries retained, an infinity of 
solutions is available. The second possible solution 
is one with a temperature maximum at x = 1/2, 
and a vanishing temperature and conductive flux 



iv 



v 



0) at X = 1. The third solution is 



one with a temperature maximum at x = 1/3, a 
minimum at x = 2/3, and again a maximum at 



X = 1, and solutions for higher eigenvalues are in 
the same vein. The segments between the temper- 
ature minimum and maximum are each still the 
inverse of the regularized /3-function, with x ^ -x 
for the segments with decreasing temperature. 

This series of solutions with an infinite num- 
ber of eigenvalues is typical for the type of second 
order ordinary differential equations that is being 
considered here, and the analogous sets of solu- 
tions for the Schrodinger, Bessel, and Legendre 
equations are well known. Note that in contrast 
to the eigenfunctions of those examples, the in- 
verted regularized /3-functions are not orthogonal, 
since they are non-negative. 

With n denoting the order of the solution, the 
expression for x of the first segment (eq. [21]) be- 
comes 



^ e„(l-2(2 + 7)/7) i/2 7/3(u;A,l/2) 
^ 2 ^ 2(2 -f 7 -ha)' 



(26) 



where e„ is the value that e must take to make this 
solution possible (i.e., the eigenvalue). Applying 
the boundary condition at the end of the first seg- 
ment {x = l/n) yields the expression for e„ (erf. 
eq. Ea). 



.2(2 + 7- 



1 



n2(l-2(2 + 7)/7) 



7 



' S2(A + 1,1/2) • 
(27) 

The physical interpretation of this result will be 
discussed in 



3. Physical Interpretation 

Having developed a simple analytical expres- 
sion for the temperature profile of an elementary 
coronal strand for a whole array of heating mech- 
anisms, I will now consider the physical interpre- 
tation of this result. First I will derive the scal- 
ing laws in their better known dimensional form. 
Then I will explore the parameter space of the so- 
lutions for different heating functions, and discuss 
how loop observations may serve as a diagnostic 
for determining the heating function. The specific 
example of Joule heating will be treated in detail, 
followed by an analysis of three other coronal heat- 
ing mechanisms that are captured by the formal- 
ism of this paper. Finally I will derive a solution 
in which the assumption of constant loop cross- 
section is relaxed and show that only marginally 
affects earlier conclusions. 



3.1. Scaling Laws 

As was pointed out in the previous section, tlie 
differential equation for the temperature profile of 
a coronal strand is over-constrained by its four 
physical boundary conditions. The consequence 
is that solutions only exist for a specific value of 
each of the two parameters (e and ^), that was 
determined above. Directly therefrom the scaling 
laws are derived. Writing out the relation for e, 
equation (|23|) . yields 



PnL — Trr 



^Xo 4 + 27 + 2a ^ 



1,1/2). 

(28) 

Note that is very close to the original RTV scaling 
law for uniform heating. For 7 — 1/2 one finds 
PqL ~ T^, the original result. The only place 
where the functional form of the heating function 
plays a role is in the constant of proportionality, 
but there it will be found to be of great potential 
significance as a diagnostic. Craig et al. (1978), 
and later Kuin & Martens (1982), derived that the 
form of the scaling laws is preserved for the non- 
uniform heating considered here, and Bray et al. 
(1991, p. 297) calculated the constant of propor- 
tionality explicitly. Kano & Tsuneta (1995, 1996) 
included a non-zero conductive flux at the loop 
apex in order to model looptop heating - as one 
would expect in flares - therefrom deriving an up- 
per and lower limit to the constant of proportion- 
ality in the first scaling law. 

The expression for A is given in equation (jlSp . 
Expressing this relation in terms of the parameters 
in the heating function yields 



A + 1 



11/2-^7 
2(2 + 7 -ha)' 



(29) 



Thus for 7=1/2, A+l=6/(5+2a), which is a de- 
creasing function of a. The first scaling law, equa- 
tion (|^5|) . demonstrates that for a loop with a 
given length and maximum temperature the loop 
pressure, Pq is proportional to B{X + 1, l/2)/(4 -|- 
27 -I- 2a) , and Figure [2] shows the degree of over- 
pressure for loops with heating concentrated at 
the footpoints (a <0), and v. v. the under-pressure 
of strands heated at the apex. Note that over- 
pressure with respect to the uniform heating is 
therefore consistent with the scaling laws, and in- 
deed quantified by them. Over-pressure in loops 
with footpoint heating is a feature that has been 



reported frequently from hydrodynamical simula- 
tions of coronal loops (e.g., Serio et al. 1981; As- 
chwanden et al. 2001; Miiller, Hansteen, & Peter 
2003; Miiller, Peter, & Hansteen 2004; Miiller et 
al. 2005). 



2 - 




a 



Fig. 2. — Over-pressure with respect to the RTV 
scaling law for uniform loop heating as a func- 
tion of the parameter a in the non- uniform heating 
function of this paper. Eh = HPqT". 

I further note that the measurement of over- 
pressure in elementary strands (as defined in the 
introduction), is, in principle, a diagnostic for the 
form of the heating profile, at least for negative 
a's, since it directly leads to an estimate of a. 
However, when a bundle of elementary strands fills 
an observed loop with an unknown filling factor, 
the emission measure of the loop cannot, by it- 
self, lead to a correct estimate for the density, and 
hence the pressure (e.g., Martens, van den Oord, 
& Hoyng 1985, Porter & Khmchuk 1995). Spec- 
tral information from density sensitive line ratios 
may be the solution here, but then there still may 
be confusion from the linc-of-sight addition of sev- 
eral strands. Below I will show that the expan- 
sion of the loops from footpoint to top (cfr. De- 
Forest, 2007) does not infiuence the over-pressure 
in a significant way, and hence strand pressure 
is a good diagnostic for the heating function. A 
word of caution must be added here though, and 
that is that this result only applies to quasi-static 
strands. Warren, Winebarger, & Hamilton (2002) 
and Winebarger, Warren, & Seaton (2003) have 
demonstrated that strands that are cooling down, 
e.g., in the aftermath of a sudden heating event 
like a flare, tend to have over-pressures with re- 
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spect to the scaling law for uniform heating as 
well, and by a far greater degree than footpoint 
heated loops 

The second scaling law is derived by writing 
down the relation for ^, equation (|11|) . explicitly. 
In terms of the original parameters for the heating 
function the result is 



max 
2-/5 



Xo^o 



7/2 + a 
3/2-7- 



(30) 



For uniform heating, a=(3—0, one finds back the 
RTV second scahng law, H = Eh ~ P'^l^jU'l^, 
after elimination of T^ax with the use of the first 
scaling law. A more intuitive way of expressing the 
second scaling law above is in terms of the heating 
rate at the top of the strand. 



E, 



top 



max 



,0 _ Fo'xo(7/2 + «) 



TlVx{3/2 - 7) 



(31) 



The two scaling laws completely define the quasi- 
static hydrodynamic solution for the strand, i.e., 
one can algebraically derive the loop pressure and 
temperature profile from the heating function - 
the apex heating rate plus parameters a and /? - 
and the loop length. 

Finally, instead of normalizing the heating func- 
tion by the apex heating, one might want to ex- 
press this function in terms of the Poynting flux 
F through the footpoints. This is accomplished as 
follows 

F = f Ehis)ds = HLP^T^^^^ f edx. (32) 
Jo Jo 

Replacing dx using the differential form of equa- 
tion (fT7|) yields an integral that can be trans- 
formed into a P function in the same manner as in 
Section [2 The final result is 



F = [ Ehis)ds = HLPI^T, 
Jo 



„ B{a + 1,1/2) 
^S(A + l,l/2)' 
(33) 
where the parameter a is defined as 



cr = A + 



a + 7 + 2 



(34) 



The requirement a > -1 translates into a > - 
(2-I-7), which is the same requirement as for A. 

Equation ([55]) above can be used to eliminate 
the constant H in the heating function in favor 



of the, in many cases, physically more meaning- 
ful Poynting flux entering from the chromosphere. 
Kano & Tsuneta (1995) have derived such an ex- 
pression without solving for the temperature pro- 
file. After eliminating H the expressions for the 
hydrodynamical variables Tmax and Pq in terms 
of the magnetic field parameters loop-length and 
Poynting flux, as well as heating function param- 
eters a and f3 become 



T 

J- rr. 



4(2 + 7)^ 12/7 



^«o^ \7/2 + a)BxB, 



(35) 



and 



Po 



^(ll+27)/142^(27-3)/14 
,(2+7)/7,l/2 



X (36) 



^0 



Xo 



2(2 + 7 + a) ^'^+^^/\ (3-27)i3A un+27)/i4 
(3-27)1/25^ \7 + 2a)Bj 



IS 



where the meaning of the shorthand Ba and B, 
obvious (cfr. eq. [33]). 

Note that only the expressions between square 
brackets depend on the details of the heating dis- 
tribution, while the parameter f3 from the heat- 
ing function does not appear at all. One inter- 
esting inference from this result is that the loop 
pressure is almost proportional to the Poynting 
flux entering into and dissipated in the loop, but 
only very weakly dependent on the loop length (~ 
L~^"). The maximum temperature of the loop is 
not strongly dependent on either: for a loop with a 
given length an increase of incoming flux of an or- 
der of magnitude only results in slightly less than 
a doubling of the temperature. Hence T^ax is a 
relatively insensitive diagnostic for the incoming 
flux, unlike the loop pressure. The fact that there 
is a large spread in observed loop apex tempera- 
tures, from about 1 to 5 MK, then implies that 
there is an enormous range in incoming Poynting 
flux for solar coronal loops. 

3.2. Temperature Profiles as a Heating Di- 
agnostic 

One might expect that vastly different heating 
functions (e.g., loop apex heating versus footpoint 
heating) would give rise to temperature profiles 
that are sufficiently different from each other to 
allow for the determination of the heating profile 
from a series of temperature measurements along 
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the axis of an elementary strand. This turns out 
not to be the case as is demonstrated in Figured 
The left hand panel of this figure shows the run of 
the dimensionless temperature for ten values of a, 
defined by 



a{n) =-2.0 + n/2; n = 0,l, ,9. 



(37) 



It is immediately obvious that the differences be- 
tween the temperature profiles are very small, and 
that the inversion (i.e., reproducing the heating 
function from the temperature profile) would be 
very hard to do with real data. The flattest pro- 
files (relatively) are those for heating concentrated 
at the footpoints, but their flatness is never equal 
to that derived in the numerical solutions of As- 
chwanden et al. (2001). An investigation into the 
reason for this discrepancy will be reported in a 
future paper. 

The right hand panel of Figure [3] shows the 
heating functions for the same values of a as in the 
left hand panel and the equation above. The heat- 
ing functions for a > are normalized by their 
apex value, while for a < the heating is normal- 
ized at a footpoint temperature to = 0.01, i.e., 1% 
of the maximum temperature at the apex. The 
panel shows that the heating becomes more and 
more concentrated near the footpoints as a be- 
comes more negative. There is no analytical so- 
lution for a— -2.5 or less, and the stability of the 
solutions with a < -2 is doubtful. The results 
of a full numerical stability analysis are currently 
under investigation and will also be reported in 
a future paper. The solution for uniform heating 
is stable and within a few percent of the numeri- 
cal solution for the full radiative loss function and 
including gravity, as reported in Martens et al. 
(2000). 

Hence, the conclusion of this section is that, in 
contrast to the scaling laws, the temperature pro- 
files of elementary strands are of limited value for 
the determination of the heating profile. The ex- 
ception to that might be the case where the heat- 
ing is strongly concentrated near the top of the 
loop. This situation is captured in a basic sense in 
the solutions presented by Kano & Tsuneta (1995, 
1996) who chose a non-zero conductive flux at the 
top of the loop in order to simulate a large heat 
influx from the top, a situation that may apply 
in the case of flaring loops. Priest et al. (1998, 
2000) present approximate solutions for the tem- 



perature profile for different heating profiles and 
claim therefrom that apex heating can be dis- 
cerned from uniform and footpoint heating from 
the temperature profile. Mackay et al. (2000) dis- 
pute this by demonstrating that the errors in the 
measurements are too large to make a meaningful 
distinction. 

Here I have simulated heating concentrated 
around the loop apex by using a large value for 
the heating parameter a. The dashed line in Fig- 
ure [3] represents the solution for the case a — 25. 
It clearly exhibits a peak in the heating profile 
around the temperature maximum, but the tem- 
perature profile is only slightly different from the 
other solutions presented in the figure. I conclude 
therefrom that in principle one might indeed dis- 
tinguish between apex heating from uniform and 
footpoint heating using the temperature profile, 
but that in practice the measurement errors will 
likely interfere, which agrees with both Priest et 
al. and Mackay et al. However, heating concen- 
trated near the loop apex may apply to flaring 
loops, but, as will be shown below, typical heating 
mechanisms for non- flaring Active Region loops do 
not have large values of the heating parameter a. 

3.3. Joule Heating 

As a specific example of a physical heating 
mechanism that is captured by the formalism de- 
veloped above I will now consider Ohmic dissipa- 
tion of DC coronal currents. Joule heating is de- 
scribed by 



Eh = Vsj 



(38) 



where j is the coronal current density, and rjs the 
Spitzer resistivity. (The subscripts "S" and "0" 
are used here to avoid confusion with the vari- 
able r] used throughout this paper). The Spitzer 
resistivity depends on the local hydrodynamical 
variables as (Spitzer 1962; Allen 1973, p. 50) 



rjs = 77oln(A)T-3/2 



(39) 



with A denoting the plasma parameter, and 770 a 
constant. Although A (the number of particles 
in a Debye sphere) depends on the local density 
and temperature (e.g., Sturrock 1994, p. 11), I will 
ignore that dependence since it matters little after 
taking the natural logarithm. For typical coronal 
conditions, T = 2 MK and density n = 3 x 10*^ 
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Fig. 3. — Left: Dimcnsionless temperature solutions for the non-uniform heating function Eh = HPqT", 
with parameters from equation (P7|) . The dashed hne corresponds to a = 25. The fowest value of a 
corresponds to the flattest curve (i.e., the upper one at t=0.5), and the temperature profiles become less 
flat as a increases. Right: The dimensionless heating functions corresponding to the temperature profiles 
on the left. For negative values of a the functions are normalized at t=0.01, and for positive values of a at 
the apex. The lowest value of a corresponds to the largest ratio of footpoint to apex heating, with the ratio 
monotonically decreasing for increasing value of a. The dashed line represent the heating profile for a = 25. 



cm •^, I find /n(A) = 16.7. rfo is equal to 9.0 x 10 ^ 
esu (Allen 1973, p. 50). 

In terms of the generalized heating function 
used above one finds a = -1.5, /3 = 0, and H = 
?7oJ^'"'(A). In a force-free coronal flux tube with 
constant cross-section the current density j is in- 
deed constant along the tube, so H is nearly con- 
stant. The solution for the temperature profile and 
the heating function is shown by one of the curves 
in Figure [3] (the second flattest temperature pro- 
file, and the second-most peaked at the footpoints 
heating function) . The amount of over-pressure in 
a Joule heated coronal strand is a factor of about 
1.5 according to Figure^ That results in an excess 
brightness of a factor 2.2 compared to a uniformly 
heated strand. 

The heating rate in coronal loops that is in- 
ferred from the observed radiation and the scaling 
laws is of the order of 10~^~^ ergcm~'^s~^. Us- 
ing the value of r/o given above and a temperature 
of 2 MK, I find that a current density of 7.1 x 10^ 
statamp cm~^ is required. That corresponds to an 
electron drift velocity of about 50 km/sec, which 
is indeed much less than the coronal electron ther- 
mal velocity, the condition for Spitzer resistivity to 
apply when the ion temperature equals the elec- 
tron temperature. Note that down the legs of the 
loop the ratio of the drift velocity to the electron 



thermal velocity decreases further, since the drift 
velocity scales as T in a strand with constant pres- 
sure and current density, while the electron ther- 
mal velocity scales only as VT. This implies that 
Spitzer heating will reach its limit first at the loop 
apex, so that anomalous resistivity, leading to fur- 
ther enhanced heating, will start at that location 
and expand throughout the corona before it pen- 
etrates the chromosphere. 

It needs to be pointed out that the current den- 
sity derived above is much higher than the cur- 
rent density inferred from vector magnetograph 
observations. Pevtsov et al. (1994, 1995, and pri- 
vate communications) report an average unsigned 
force- free a per pixel of the order of 10~^ cm~^ 
in active regions. This is an order of magnitude 
larger than their average a in active regions, be- 
cause a's of both signs are present. Each ac- 
tive region appears to be subdivided in a hand- 
ful of patches with same sign a (Pevtsov, private 
communication). Since the force- free a is con- 
stant along a field line, it has the same value in 
the corona. With an extrapolated coronal field 
strength of the order of 100 Gauss one finds from 
Ampere's law j = 240 statamp cm~^, more than 
four orders of magnitude shy of the value required 
for Ohmic heating derived above. 

However, the measurements of Pevtsov et al. 
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at the Haleakela Stokes Polarimeter are samples 
with a pinhole of approximately a six arcsecond 
diameter, and the currents derived therefrom are 
totals over the whole pixel, not really current den- 
sities. Hence it is possible, in principle, that there 
are sheets with very high current densities in each 
pixel, with nearly current-free regions in between, 
that produce the required heating rate, although 
at the obvious expense of a reduced filling factor. 
It is easily verified that in a 100 Gauss field a sin- 
gle sheet with a thickness of the order of 340 m 
and a current density of 7.1 x 10^ statamp cm"^, 
as derived above, would already produce a cur- 
rent density of twice the 240 statamp cm~^ per 
pixel implied from the Pevtsov et al. results. The 
filling factor can be enhanced by introducing a 
larger number of sheets with counter-flowing cur- 
rents. The work of van Ballegooijen (1986) on cur- 
rent filamentation resulting from regular footpoint 
photospheric motions, and the work by Longcope 
(1996) on current generation on magnetic separa- 
tors, both strongly suggest that coronal currents 
will be highly non-uniform. Hence coronal heating 
by Ohmic dissipation in a number of unresolved 
high current density sheets inside coronal loops 
is a candidate heating mechanism that cannot be 
dismissed out of hand. 

3.4. Anomalous Resistivity Phase-mixing, 
Nano-flares, and Velocity Filtration 

It was mentioned above that a number of coro- 
nal heating mechanisms proposed in the literature 
can be captured by the heating parametrization 
used in this paper, with the heating rate propor- 
tional to density (or pressure) and temperature, 
each to an arbitrary power (eq. [2]). In this sec- 
tion I will briefly describe a number of these in 
addition to regular Joule heating with Spitzcr re- 
sistivity discussed above. This discussion does not 
constitute a comprehensive review. 

First let's consider anomalous current dissipa- 
tion driven by the ion-acoustic instability, which 
occurs when the electron drift velocity surpasses 
the electron thermal velocity of the plasma. In 
that case the resistivity is anomalously enhanced 
over the Spitzer resistivity and the current dissipa- 
tion and heating rate can greatly increase. Rosner 
et al. (1978a) discuss this scenario in depth and 



cite a resistivity of the order of 

77 = 2 X IQ-^n-^/^s 



(40) 



where n,, is the electron density in cgs units. As 
shown above the current density j is constant 
along the strand and hence, using the gas law, 
the heating rate Eh = ?yj^, scales as T^/'^P~^/^, 
Therefore, in the formalism of this paper a — 1/2 
and (3 — -1/2 for ion-acoustic heating. 

Viscous dissipation of shear Alfven waves is a 
viable heating mechanism proposed by Heyvaerts 
& Priest (1983). The physical basis is that stand- 
ing Alfven waves on neighboring field lines will 
rapidly grow out of phase even for a slightly vary- 
ing Alfven speed across the magnetic field, a pro- 
cess called phase-mixing. This produces large ve- 
locity gradients that result in viscous dissipation 
of the wave velocity fields. The heating rate in the 
Heyvaerts and Priest mechanism is given by 



Eh 



^J-vPVu 



JS 



(41) 



with [i^ representing the shear viscosity, p the den- 
sity, Vwave the velocity amplitude of the Alfven 
wave, and d the typical cross-field length-scale 
over which the wave velocity varies. Heyvaerts 
and Priest demonstrate that the cross- field length- 
scale d decreases during the process of phase- 
mixing until the total wave dissipation rate bal- 
ances the incoming Poynting flux from the Alfven 
waves. 

It may be assumed that the kinetic energy den- 
sity in the standing Alfven waves is constant along 



the loop, i.e., pv^^ 



constant. There is in- 



deed some observational support for this in that 
the measured square of the non-thermal broaden- 
ing velocity in Active Region emission lines scales 
linearly with temperature (Saba & Strong 1991), 
although these data from the Flat Crystal Spec- 
trometer onboard the Solar Maximum Mission 
{SMM) have very low resolution (~ 15 arcsec) so 
that individual loops are not resolved. Combining 
the gas-law with energy equipartition in the wave 
iP^wave — constant) reproduces this scaling. Us- 
ing the constancy of d and pw^ane throughout the 
loop I find that the heating scales with the shear 
viscosity /z. Since p depends on plasma density 
and temperature as T'°l'^ (e.g., Choudhuri 1998), 
the heating rate from the phase-mixing model is 
characterized by a = 5/2 and /3 = 0. 
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Note that according to equation ((28l) and Fig- 
ure [21 for a loop strand with a given (observed) 
length and maximum temperature the gas pres- 
sure for the case of Joule heating exceeds that for 
the case of heating via phase-mixing by a factor 
two. Hence a measurement of strand pressure, e.g. 
by taking measurements of density sensitive line 
pairs at a given temperature, is a possible diag- 
nostic for discriminating between these two heat- 
ing mechanisms, once the elementary strands have 
been resolved. 

Uniform heating, i.e., a = and (3 — 0, is the 
type of heating often considered in numerical sim- 
ulations, and also in the original scaling-law pa- 
pers by Rosner et al. (1978b) and Craig et al. 
(1978) . The reason is of course that this is the sim- 
plest assumption, but also that the widely known 
theory of nano-flare heating of the solar corona 
(e.g., Parker 1988, 2005) is compatible with uni- 
form heating. Because of the high plasma /3 in 
coronal loops the origin and development of nano- 
flares is considered by Parker as a purely mag- 
netic phenomenon, i.e., not dependent in first or- 
der upon local thermodynamic quantities such as 
density and temperature. Hence it is reasonable 
to assume that the distribution of the nano-flares 
is uniform throughout loop strands. Moreover, in 
order to reproduce the observed quasistatic ap- 
pearance of Active Region loops (Antiochos et 
al. 2003), for each loop strand the time between 
nano-flares has to be shorter than thermal cool- 
ing time of the strand. And indeed, because the 
loop thermal cooling time acts as an effective fre- 
quency filter in the hydrodynamical evolution of 
loops, nano-flares with a frequency that is suffi- 
ciently smaller than the frequency associated with 
the cooling time, will appear, to first order, as 
steady heating. Of course there will be effects such 
as Doppler-shifts and line-broadening due to the 
explosive nature of the nano-flares., which can be 
used as diagnostics to distinguish between nano- 
flare heating and true steady heating (e.g., Pat- 
sourakis & Klimchuk, 2006) but the thermal struc- 
ture of the strands will be identical to first order. 

Finally, some heating mechanism are consid- 
ered non-local, in the sense that free magnetic 
energy is converted into particle kinetic energy 
at a different location from where particle kinetic 
energy is thermalized in the background plasma. 
Thus the conversion of magnetic to particle en- 



ergy takes place in a different location than the 
heating of the background plasma, and hence the 
term non-local, although the heating itself still oc- 
curs in-situ. One example of this is the heating of 
flaring loops resulting from the impact of electron 
beams that are often assumed to have been accel- 
erated in a current sheet on top of these loops. A 
second example, of potential relevance to the heat- 
ing of quasi-steady Active Region coronal loops, is 
the filtration of the high energy non-thermal tails 
in the electron distribution from the transition 
region into the overlying corona (Scudder 1992a, 
1992b, 1994; Vihas, Wong, & Khmas 2000). We 
then have a situation in which high energy elec- 
tron beams enter coronal loops through their foot- 
points, emanating, paradoxically, from the much 
cooler transition region and chromosphere. Scud- 
der points out that this effect is possible because 
of the much higher particle density in the chro- 
mosphere and transition region. Since the mean- 
free-path of the electrons in the filtration beams 
scales with the fourth power of the beam electron 
velocity, it only takes a relatively small ratio of 
beam velocity to thermal electron velocity in the 
corona for the mean-free-path of the beam elec- 
trons to surpass the loop length and hence for the 
filtration beams to travel through the loop strands 
relatively unattenuated. In that case the number 
of collisions between beam and background plasma 
particles at any location in the loop is proportional 
to the density alone, and therefore the heating rate 
is as well. In the formalism of this paper that re- 
sults in the parametrization a ~ -1 and /3 = 1. 

In closing this section I note that the coronal 
heating mechanisms discussed here (anomalous 
current dissipation, phase-mixing, nano-flares, 
and velocity filtration) are all relatively old. The 
reason is, I believe, that after it became clear that 
coronal heating through the dissipation of acous- 
tic waves was implausible (Athay & White 1979), 
solar theorists picked the low-hanging fruits first, 
i.e., the most obvious and plausible heating mecha- 
nisms were investigated first. More current papers 
on coronal heating often continue the development 
of older ideas, some of which are described here, or 
pursue more exotic possibilities. Solar physics is 
now faced with the situation of a plethora of pos- 
sible coronal heating mechanisms, which, as the 
current paper illustrates, are difficult to distin- 
guish observationally since the temperature pro- 
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files are almost identical, and the predicted loop 
pressures vary by amounts that are currently dif- 
ficult to derive from measurements of unresolved 
multi-stranded loops. However, given this mul- 
titude of theoretical possibilities, I disagree that 
our inability to pinpoint the actual coronal heat- 
ing mechanism forces the conclusion that there is 
no in-situ coronal heating at all (Aschwanden et 
al. 2007)! 

3.5. Loops with Variable Cross-Section 

As was pointed out in the introduction, a re- 
cent paper by DeForest (2007) has made a pow- 
erful argument that the elementary strands of at 
least EUV emitting loops expand with distance 
along the loop axis, as one would expect based 
on simple potential or force-free magnetic field ex- 
trapolations. The formalism developed in this pa- 
per can be used to capture loop expansion with 
height in a straightforward way. For a loop with 
variable cross-section the energy equation (eq. [l]) 
takes the form 



W)i''''^^""'i:' 



where A{s) represents the variable cross-section as 
a function of distance along the axis of the elemen- 
tary strand. To make this equation amenable to 
analytical treatment Following Bray et al. (1991, 
p. 294) I assume a power-law relation between 
cross-section and temperature. 



A{s) = Ao{T{s)/T^axY 



(43) 



where A{Q) is the cross-section at the apex of the 
strand, and ^ is a parameter. Note that for a 
strand that expands with height 5 is positive, and 
that constant cross-section corresponds to 5 — 0. 
Using the same technique as in § [2]to reduce the 
energy equation to its dimensionless form I find 
the remarkable result that the energy equation in 
terms of the variable rj (eq. [5] ) , with parameters 
V and /i remains unchanged. The expression for A 
in terms of v and /i, equation (jlSp . as well as the 
boundary conditions, remain unchanged as well, 
while the definition of the parameters v and /i be- 
comes 



a + 5 

7/2 + (5' 



(44) 



A* = 



^ - (2 + 7) 
5 + 712 



(45) 



The main difference with previous results is in the 
definition of the dimensionless temperature, 



t 



rj'' 



(46) 



which results in only a small correction in the ex- 
pression for the temperature profile of the strand: 

T{S) = T„,aAPn\s/L- A + 1, 1/2)] (i+2./7/(2+,+„, . 

(47) 
The other differences are in the expressions for ^, 
and e. After some straightforward algebra, the 
explicit expressions for the scaling laws of tapered 
fiux tubes become 



E, 



top 



Pixo{7/2 + a) 



PqL — T„ 
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J- max V'J/ ^ 
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(48) 



4 + 27 + 2a 



1 + 25)' 
S(A + 1,1/2). (49) 



Without solving for the temperature profile Bray 
et al. (1991, p. 297) have derived the same results 
for the scaling laws of tapered loops. 

The difference of the expressions above with the 
scaling laws for constant cross-section amounts to 
a correction factor 1/(1-1-2(5/7) in the power-law 
relation between t and rj for the temperature pro- 
file (eq. [47]), and a minor correction in the second 
scaling law. However, in the first scaling law the 
correction can be significant. A value oi 5 — 1, 
which as will be shown below is not extravagant, 
leads to an over-pressure of a factor -^3 ~ 1.7, 
which is of the same order as the over-pressure 
that can be generated by footpoint heating in com- 
parison with uniform heating. This complicates 
matters, but I will show below that potentially 
the integrated brightness of "moss" can be used 
to determine the amount of loop tapering inde- 
pendently. 

Figure |4] illustrates the effect of flux tube ex- 
pansion on the solutions, for the case of current 
heating with Spitzer resistivity, discussed above 
(q:=-3/2). The left panel shows the temperature 
profile for both constant cross-section {5—0) and 
a loop expanding in radius as the square root of 
temperature ((5=1, dashed curve), and the right 
panel shows the loop radius as a function of dis- 
tance along the axis from bottom to top. The as- 
sumed temperature dependence of the expansion 
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Fig. 4. — Left: Dimcnsionless temperature solutions for an Ohmically heated (Eh ~ T~^'^) elementary 
strand with constant cross-section ((5=0, thin line), and an increase in strand radius proportional to vT 
{6—1, dashed line). Right: Dimensionless loop radius as a function of distance along the loop axis for the 
expanding loop. 



rate produces roughly a factor ten ratio in loop 
radius between the apex and the chromospheric 
footpoint of the strand, in agreement with the ex- 
pansion factors that follow from field extrapola- 
tions (see DeForest, 2007). 

These results show that here is only a rela- 
tively small difference in the temperature profile 
and pressure between strands with uniform and 
expanding cross-sections. Therefore the ratio in 
loop integrated emission measure between a sec- 
tion of the loop near the top and one near the 
footpoints constitutes a good diagnostic for the 
expansion factor, since the emission scales linearly 
with the cross-sectional area. At the same time, 
for well resolved strands the emission per pixel will 
barely be influenced by the expansion factor. 

Note that strand expansion may have a signif- 
icant influence on the emission from moss, foot- 
point EUV emission from hot loops (Peres, Reale, 
& Golub 1994; Schrijver et al. 1999; Berger et al. 
1999; Fletcher & DePontieu 1999). It was pointed 
out by Martens et al. (2000) that the moss emis- 
sion can be used as a diagnostic for the strand 
footpoint constriction. In the example of Figure |4] 
a hot strand with a maximum temperature of 3.5 
MK, will have its cross-section area constrained by 
a factor 0.29 at 1 MK, the temperature that dom- 
inates the emission in the TRACE 171 A channel, 
which shows the moss very well (e.g., see Fig. [T]). 
Since according to Figure[4]there is little change in 
the temperature gradient at that point, the emis- 



sion measure of the moss is also reduced by a factor 
of approximately 0.3 when compared to a loop of 
constant cross-section. 

Fletcher & DePontieu (1999) and Martens et 
al. (2000), both using TRACE data, as well as 
recently Warren et al. (2007, private commu- 
nication) and Winebarger, Warren, & Falconer 
(2008), using the Hinode EUV Imaging Spectrom- 
eter (EIS) find a filling factor for the moss of the 
order of 10%, less than the result given above. 
As discussed in Martens et al. the filling factor 
of the moss is the product of the filling factor in 
the coronal portion of the loop times the magnetic 
constriction near the loop footpoints. The calcu- 
lation above only considers the magnetic constric- 
tion factor, and combined with a coronal filling 
factor also of the order of 0.3 would lead to the 
observed result of about 10%. Independent obser- 
vations of the filling factor in the coronal portion 
of the loops are needed to separate the two com- 
ponents of the filling factor. For the purpose of 
estimating the constriction factor (5, as discussed 
above, only the ratio between the loop apex emis- 
sion measure and that of the moss needs to be 
known. 



Solutions with Temperature Rever- 
sals: Physical Interpretation 



3.6. 



In § 12.21 the existence of solutions with tem- 
perature reversals has been demonstrated. The 
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question arises whether such solutions can occur in 
physical reality, and I will demonstrate here that 
in certain circumstances the answer may be posi- 
tive. 

In the simplest solution with a temperature re- 
versal {n=2) one would find a condensation with a 
chromospheric temperature at the top of the loop. 
Such solutions have been investigated numerically 
by Karpen et al. (2001, 2003, 2005, 2006), and 
these authors found that they are unstable in the 
usual semi-circular loops because of gravity: the 
cool "blob" will be pulled down, in a similar fash- 
ion as in the Rayleigh- Taylor instability. It is rea- 
sonable to expect the same effect for higher order 
solutions with multiple cool blobs. 

However, in order to investigate plausible ge- 
ometries for solar prominences, Karpen et al. fur- 
ther simulate strands that are stretched out in 
the horizontal direction (i.e., loops that are mostly 
horizontal, except near both footpoints), and also 
loops that have a dip (a U-shaped section) in 
the middle, a widely accepted physical model for 
sustaining prominence material in the corona far 
above its scale-height. They find that such situa- 
tions can be thermally stable in some cases, or are 
maintained for relatively long times before dissi- 
pating in others. I suggest that these numerical 
solutions are a manifestation of the temperature 
reversal solutions derived in this paper. 

The following thought experiment further ver- 
ifies the physical viability of the temperature re- 
versal solutions. The solution with a uniformly in- 
creasing temperature stretches from one footpoint 
to the top of the loop. In the other half of the loop 
the solution must be the same, i.e., same T^ax, 
loop half-length L, and pressure Pq, but mirrored 
in the space coordinate of course. That situation is 
physically identical to the n=2 solution for a loop 
with the same maximum temperature but twice 
the loop length. In particular the loop pressure 
should be the same. 

By writing out the eigenvalue solution for e„ 
(eq. [27j ) in terms of the original physical param- 
eters one finds the scaling law(s) for loops with 
temperature reversals, 



PoL/n = T, 



^r^iV^tM^i?(A+i,i/2). 



Xo 



A + 2-f + 2a 



(50) 
This shows that the n=2 solution for the loop with 



twice the loop-length as the original one, and with 
the same temperature and heating profile, indeed 
produces the same loop pressure, as the thought 
experiment requires. This scaling law also demon- 
strates that for the same loop with the same heat- 
ing profile and maximum temperature, the loop 
pressure in the case of temperature reversal pro- 
files is an n-multiple of the pressure of that in the 
loop without condensations. 

I conclude from this section that it is plausi- 
ble that the formalism of this paper captures not 
only the temperature profile of ordinary hot coro- 
nal loops, but also, in a crude sense, that of solar 
prominences carrying condensations. 

4. Discussion 

Based upon the observations that emission 
from the solar corona is dominated by Active Re- 
gion soft X-ray loops, and that these loops are 
mostly stable for times longer than their ther- 
mal timescales (e.g., Antiochos et al. 2003), and 
coupled with the likelihood that observed coronal 
loops consist of a multitude of thermally isolated 
strands, I have set out to derive analytical mod- 
els for the thermal structure of quasi-static and 
isobaric coronal strands. I have derived analytical 
solutions for the temperature structure and scal- 
ing laws in coronal strands in a formalism for the 
coronal heating profile that captures loop foot- 
point heating, uniform heating, as well as loop 
apex heating. Solutions with both uniformly in- 
creasing temperature and with temperature rever- 
sals have been derived, and it is suggested, from 
a comparison with numerically derived results in 
Karpen et al. (2006, and references therein) that 
the latter solutions apply to the temperature pro- 
file in prominences. 

The analytical solutions produce strand tem- 
perature profiles that vary very little for enormous 
differences in the distribution of the heating over 
the loop, thus making the temperature profile a 
rather poor diagnostic for the heating functions. 
On the other hand, I demonstrated that the con- 
stant of proportionality in the RTV scaling law 
(PqL ~ T'^aa;) dcpeuds ou the distribution of the 
heating over the strand, with footpoint heating 
producing over-pressure in comparison to uniform 
heating. Hence, measurements of pressure or den- 
sity, for instance from density sensitive line ratios. 
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in combination with temperature measurements 
may be very useful in teasing out the distribution 
of the heating over observed loops. 

The analytical solutions have been derived for 
strands with uniform cross-section, but a gener- 
alization to a family of solutions for non-uniform 
cross-sections has been demonstrated as well, and 
it turns out that these solutions only differ only 
by a small correction factor from the solutions 
for constant strand cross-section. Through a sim- 
ple example it is demonstrated that variation in 
strand-cross section can account for the deficit in 
observed "moss" emission in comparison with the 
models. As has already been pointed out by De- 
Forest (2007) the same effect can account for the 
observed difference between scale-height tempera- 
ture and observed temperature in EUV loops. 

I have further shown that the expression used 
for the coronal loop heating function is quite ver- 
satile, by taking a quintet of well known leading 
candidates for the coronal heating mechanism and 
explicitly expressing them in terms of the heating 
formalism. However, I should point out here that 
certainly not all proposed coronal heating mecha- 
nisms can be expressed in the form Eh ~ PqT". 
For example, a wave flux entering a strand from 
the chromosphere, and dissipating on a length- 
scale £ that is short compared to the loop-length, 
will lead to a heating profile of a form that is ex- 
plicitly dependent on the distance from the foot- 
point s, e.g.. Eh ~ exp(— s/£), and that there- 
fore cannot be expressed in terms of local thermo- 
dynamic variables alone. Such an expression has 
been used in numerical simulations by a number of 
authors, for instance Serio et al. (1981), Aschwan- 
den et al. (2001), and Miiller et al. (2005). Their 
numerical solutions exhibit a variety of features 
that are very similar to those of the analytical so- 
lutions derived here: there is over-pressure in loops 
with footpoint heating, and when the heating is 
concentrated too strongly near the footpoints the 
quasi-static solution ceases to exist (Miiller et al. 
2005), just as I find for a < -2.5. 

The temperature profiles presented here are ex- 
pressed in terms of the inverse of the incomplete 
regularized /3-function. This type of solution has 
first been derived by Martens (1981) for applica- 
tion in stellar coronae, and has been applied to 
solar coronal loops by Kuin & Martens (1982) for 
the case of uniform heating. The solution has been 



used extensively by for the purpose of solar and 
stellar loop modeling by Fisher & Hawley (1990) 
and Hawley & Fisher (1992, 1994). To my knowl- 
edge the analytical solutions for the temperature 
profile expressed as incomplete /9-functions in the 
current paper, are the first for non-uniform heat- 
ing profiles. As pointed out above, Craig et al. 
(1978) were the first to show that the form of the 
first scaling law is preserved for non-uniform heat- 
ing expressed in terms of pressure and density to 
a given power (eq. [2]). Bray et al. (1991) first cal- 
culated the heating dependent constant of propor- 
tionality. Kano & Tsuneta (1995, 1996) added a 
non-zero incoming thermal conductive flux at the 
loop apex to model the apparent loop top heat- 
ing in flaring loops and derived expressions for 
the scaling laws in terms of inequalities as well 
as a scaling law for loops heated completely at the 
apex. 

As mentioned in the introduction, having an- 
alytical solutions for the temperature profile and 
pressure in loop strands for a number of leading 
candidate coronal heating mechanisms makes it a 
lot simpler to calculate DEMs per pixel in forward 
modeling of coronal loops consisting of a multitude 
of unresolved strands. This kind of capability will 
be most useful for analyzing coronal loops with the 
AIA instrument onboard SDO that will observe 
coronal loops through six narrowband EUV band- 
passes, in combination with XRT on Hinode that 
has nine broadband filters for observing soft X-ray 
emission. However, before this capability can be 
used one has to investigate the existence, stability, 
and accuracy of the analytical solutions derived 
here. In particular one has to analyze the effects 
of the assumptions of constant pressure, zero tem- 
perature and flux at the lower boundary, and the 
approximation of the radiative loss function by a 
single power-law. For the case of uniform heat- 
ing {a=f3—Q in the heating formalism of this pa- 
per), the existence and stability of the solutions 
has been confirmed many times over by numerical 
simulations. Martens et al. (2000) also considered 
the accuracy of this solution by comparing the nu- 
merical and analytical solution, and they found a 
temperature difference of a few percent at most at 
every location in the loop for loop heights smaller 
than the pressure scale height for the apex tem- 
perature. 

The numerical analysis of the solutions for non- 
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uniform heating is in progress (Martens, Winter, 
& Munetsi-Mugomba 2007) and has yielded so far 
that for a > — 2 stable static solutions exist, and 
that are very accurate, as long as the loop height 
is less than the apex temperature pressure scale 
height. For heating that is strongly concentrated 
near the footpoints we found cyclic solutions very 
similar to those obtained by Miiller et al. (2005) 
(and anticipated by Kuin & Martens, 1982) for 
very small values of the dissipation length £ (see 
above) compared to the loop length. These results 
will be demonstrated and analyzed in much more 
detail in an upcoming paper (Winter, Martens, & 
Munetsi-Mugomba 2008, in preparation). I note 
that I have found the study of stability through 
an analytical linear stability analysis too cumber- 
some, and inevitably leading to a numerical analy- 
sis anyway, which is why we attacked the problem 
with a full time-dependendcnt numerical simula- 
tion from the onset. 

In conclusion then, I have generalized the an- 
alytical solutions for the scaling laws and tem- 
perature profile first published by Rosner et al. 
(1978b) from the case of uniform heating to a 
family of non-uniform heating functions, includ- 
ing footpoint and loop apex heating. These results 
lead to increased insight in the relations between 
the thermodynamic variables and heating parame- 
ters in loop strands, and will substantially reduce 
the effort required for forward modeling of pixel 
Differential Emission Measures of multi-stranded 
loops to be compared with combined data from 
XRT and AIA. 
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A. Integration for Non-zero Chromospheric Temperature 

The dimensionless integral expression for the temperature profile in equation (|f7|) in § 12.11 is 

[e(^+l)/2]i/2 r z^dz 



M 7o [i-zY'^' 



(Al) 



with the parameters defined in that same section. When the temperature at the chromospheric boundary 
at a; = is not ignored, we a have lower boundary condition uo defined by uo = [ ^fe'-°m j2+a+7 which is a 
very small number except for a [ —2.5. 



Equation (jAip becomes 

[e(Ai+l)/2]i/2 r" z^dz p z^dz 



or, after integration 

^ e(l - 2(2 + 7)/7) 1/2 7[/3(^.; A + 1, 1/2) - /^K; A + 1, 1/2)] 
"^ ^ 2 ^ 2(2 + 7 + a) 



(A2) 



(A3) 



For the same value of e the extra term would merely amount to a small shift to the left along the x-axis for 
the solution. However, because the boundary condition u = 1 at a; = 1 also has to be met, the value of e - 
and thereby the constant in the first scaling law - will change. 
Applying the boundary conditions at the top of the loop yields 

^ 6(1 - 2(2 + 7)/7) .i/2 1[B{\ + 1, 1/2) - /3(^o; A + 1, 1/2)] 

^ 2 ^ 2(2 + 7 + a) ' ^ ' 

where again B is the complete j3 function. The value of the term expressed by the incomplete /3-function is 
small compared to that of the complete /3-function, B(A+l,l/2) for A > -1, i.e., a > -2.5, as will be shown 
presently. 

For small uq the denominator in the integrand term of the incomplete /3-functioncan be set to unity, i.e.. 

Restoring the original parameters t^, a, and 7 (see § 12. ip yields 

uo^+i (2 + a + 7)C/^+'^/^ 



A + 1 11/4 + 7/2 



(A6) 



which is of the order of 10^^ for Iq ~ 10^^, while the value of the complete /3-function is of the order 0.1 
for the allowed positive values of A. Note that for a | —2.5, Uq is no longer very small, even while ip is, but 
despite of that the value of the incomplete /^-function remains small, because the value of A in the integral 
expression of equation (|A5|) becomes very large. This has been verified numerically. 

Hence equation (jA4|) produces a correction of the order of 0.3% for the value of e, and therefrom a 
correction of a fraction of again 10~^ in the loop pressure according to equation (|28p. a correction that I 
have ignored in the body of the paper. 
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